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Problem (0) 


A program SMARAND has been designed to generate S(n) up to a preset limit N (N up to 
1000000 has been used in some applications). The program requires an input of prime 
number up to VN. Initially the program caculates S(p;**)=D(i,k) for all primes p; and all 
exponents k needed to reach the preset limit. It then proceeds to factorize consecutive 
values of n. If n is prime then S(n)=n otherwise p;" is replaced by D(i,k) whenever k> 1. The 
largest component in the resulting array is determined and is equal to S(n). Slighly different 
versions of the program has been used depending on the application. Up to n4 32000 both 
S(n) and D(i,k) were registed on files. The values of S(n) for n « 4800 were listed with the 
help of a program SN TAB and the values of S(p;*) were listed for p< 73, k < 75 with the 
help of a program SNP TAB. 


Problem (i) 

(a) No closed expression for S(n) has been found 

(b) No asymptotic expression for S(n) has been found. The behaviour of S(n) for n«32000 
has been graphically displayed using a program SN DISTR. 


Problem (ii) 
S(n)/(n-m), m£0 is equivalent to 


m=n+kS(n), k integer .... (1) 
Let us assume that m is a given prime p. From the definition of S(n) it is evident that for 


every n there exists a prime q such that q/n (or n=lq, | integer) and S(n) =jq (i integer). We 
can therefore write (1) in the form 


p=q(1+k) ... (2) 


To find solutions to (1) when m is a prime p it is therefore sufficient to chose n as a 
multiple of p which fills the condition 1+kj=1. 


In practice (as far as I have found) this means excluding from n those multiples of m(=p) 
which are divisible by primes larger than p and also cases where n-m has a different parity 
from S(n) as for example (n,m,S(n)) =(5054,19,38) is not a solution while (2527,29,38) is a 


solution. 


When m is not a prime let p be the largest prime such that p/m, i.e. m- rp. Solutions to (1) 
will then be found when n is a multipe of p for which the GCD of n and S(n) is rp. 


The above conditions are sufficient but may not be necessary. Lists of solutions are however 
easily obtained (not included) by looking for solutions to (n-m) mod S(n) = 0. 


Problem (iii) 

A number of solutions to S(x")+S(y") =S(z*) has been obtained and listed for n=3,5,7 and 
11. The program SMAR iii uses only Smarandache function values of the type S(p 
which had first been sorted in ascending order using a program SNP SORT. 


Problem (iv) l 
A program SMAR iv has been designed to find solutions to the equation S(k")'=S(n*) but 
no non-trivial solutions were found in the selected search area nx 8000. 


Problem (v) 

In a first attempt values saved on file up to 32000 were in a program SMAR v to calculate 
the sums Z(n)=1+1/S(2)+1/S(3)+ ... * 1/S(n) for n= 800, 1600, 2400, ... 32000. These 
sums were used to study the behaviour of Z(n)-T(n) for various functions T(n): 

T(n) =log(S(n)) gave a curve "parallel" to Z(n). 

T(n) =log(largest prime « n) gave a similar result. 

T(n)=1+1/2*+1/37+ ... * 1/r* gave interesting results. Supplemented with a linear term 
a "nearly straight horizontal" line was obtained. 

To see if this holds for larger values the exercise was repeated for n< 1000000. Computer 
files to store S(n) is now out of question and the generating program SMARAND was 
revised so that the partial sums Z(n) were calculated in the same program. T(n) was 
calculated in a separate program for vaious values of a. For a=0.5 and it was found that 


1+1/S(2)+1/S(3)+ .. +1/S(n) - (11/2 1/V3- ... +1/Vn) - (20k-58) 


where k= n/25000 deviates from 0 with at most 10 in the interval 1 to 1000000 (at the 
points of representations in the graph, 1000000 was divided in 40 interval of 25000). 


Problem (vi) 
Not attempted. 


Equipment 

Calculations were done on an dtk-computer with 486/33 Mhz processor. Programs were 
written in Borlands Turbo Basic. Printouts were done on an HP IIP Laser printer. Some 
graphs were done on an HP Paintjet. The run time to calculate S(n) up to 1000000 was 2 
h 50 m. The initial calculation up 32767 took 198 s. 
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'SMARAND, H. Ibstedt, 930320 


‘The Smarandache function S(n) calculated by comparing largest prime and S(P^A). The 
values of S(n) are calculated and registered in a file SN.DAT up to n = 32000. The calculation 
goes further in other applications. 


DEFLNG A-S 
CLS :T- TIMER 

DIM P(168),D(168,75),K(168),L(168) 

OPEN "PA" FOR INPUT AS #1 

FOR |=1 TO 168 :INPUT s1,P(I) :NEXT :CLOSE #1 


"This part of the program calculates S(P(I)^A) and saves the result in the array D(I,A), P(I) is the 
lth prime number. The routine uses the fact that D(I,A) < =P(I)*A for a downward search for the 
value of D(LA). This calculation goes beyond what is required to calculate S(n) up to n - 
32000. 


FOR 121 TO 42 
Az2:PzP(I) :D(,1) P 
WHILE A « 76 
C=0 :N=0 
L: 
C=C+1 
N=N+P 
IF C» =A THEN D(L,A) =N :GOTO LWEND 
PP=P*P 
L1: 
IF N-PP*INT(N/PP) =0 THEN C=C+1 :PP=PP*P :GOTO L1 
IF C» =A THEN D(A) =N :GOTO LWEND :ELSE L 
LWEND: 
INCR A 
WEND 
NEXT 


‘The array D(I,A) is stored in a file SNP.DAT for future use. 


OPEN 'SNP.DAT" FOR OUTPUT AS #2 
FOR 1«1 TO 42 :FOR J=1 TO 75 
PRINT #2,P(I),J,D(1,J) 

NEXT :NEXT :CLOSE #2 


This part of the program calculates S(N). It calls on the subroutine NFACT to express N in 
prime factor form. Factors P(I)^A with A» 1 are replaced by D(I,A) and placed in array L() 
together with the factors P(I) of multiplicity 1. S(N) is then the largest component of L(). S(N) is 
stored in a file SN.DAT. 


N=1 

OPEN 'SN.DAT" FOR APPEND AS #3 
WRITE #3,1 

WHILE N «32000 

INCR N :print n 

'Factorize N. 
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GOSUB NFACT 
IF K(0) 20 THEN S - P(0) :GOTO LWR 


‘Construct L(). 


FOR 1-1 TO 168 :L(1) =0 :NEXT 
C=0 

FORI-1TOM 

INCR C 

IF K(l) -1 THEN L(C) - P() 

IF K(I)>1 THEN L(C) = D(I,K(1)) 
NEXT 


Find the largest value of L() and hence S(N). 


S-0 
FORI-1TOC 
IF L(I)» S THEN S -L (I) 
NEXT 

LWR: 
WRITE #3,S 
WEND 
CLOSE #3 
T - TIMER-T :PRINT T 
END 


‘Subroutine for factorization of N. 


NFACT: 
FOR 1-0 TO 168 :K(I) 20 :NEXT :P(0) =0 
N1=N :M=0 
FOR |=1 TO 168 
LA: 
IF N1-P(I)*INT(N1/P(I)) =0 THEN K(I) -K(I) +1 :M=1:N1 zN1/P(I) :GOTO LA 
IF N121 THEN 1168 
NEXT 
IF N1>1 THEN P(0) Z N1 :K(0) =1 
RETURN 
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